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Abstract 

We present results for the three-loop top quark contributions to the p parameter 
in the limit of large top quark mass. The simultaneous dependence on the mass of 
the Higgs boson Mh and the mass of the top quark mt is obtained from expansions 
in the range of Mh around m t and in the limit Mh 3> m t . In combination with 
the previous result for Mh = the dependence of the p parameter on the leading 
Yukawa contributions, i.e. on mt and Mh, is well under control for all mass values 
of practical importance. The effects lead to a shift in the W mass in the order of 5 
MeV and are relevant for precision measurements at TESLA. 



1 Introduction 

The standard model predicts a strong influence of virtual heavy top quarks in the low 
energy limit on four-fermion processes through their virtual contributions to the W and 
Z boson self energies pQ. This effect allows the indirect determination of the mass of the 
Higgs boson and is generally essential for electroweak precision tests. 

The bulk of the heavy top quark corrections to the W and Z self energies can be 
collected in a deviation of the so called p parameter from the tree level result Ptree = !• 
The p parameter is usually denned by the ratio of the neutral and charged current coupling 
constants at zero momentum transfer pQ: 

■WO) m 
p = ^(o)' (1) 

where Jcc(O) is given by the Fermi coupling constant Gf determined from the /i-decay 
rate whereas Jjvc(O) is measured by neutrino scattering on electrons or hadrons. 

In leading order contributions from the Yukawa coupling of the top quark to p are 
parameterized in terms of Ap. This contributions stem from the transversal parts of the 
self energies of the W- and Z-boson: 

i . i <-DT(Q) (2) 
p i-a p ~c 2 m 2 z - np(o) ' 1 ' 
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where c 2 = M^/M| represents the cosine of the weak mixing angle, defined in the standard 
on-shell scheme. Corrections from vertex and box diagrams always involve extra powers 
of the weak coupling constant g 2 weak and are thus suppressed by powers of (^woak/^Yukawa) 2 ' 
i.e. (M w /m t ) 2 . 

In this context U^ w (0) and Ilf z (0) can be taken as the bare one-particle irreducible 
amplitudes. The divergencies are absorbed by the Higgs mass, top mass, and W boson mass 
renormalization, i.e. by expressing the bare m t , Mh, and M\y in terms of the renormalized 
quantities. In the following we will first express the final answer in terms of the MS top 
mass mt, the on-shell Higgs mass M#, and the low energy Fermi constant G F obtained 
from /x-decay The latter is introduced to absorb the vector boson mass M w . In a last 
step the transformation from the MS top quark mass m t to the on-shell top quark mass 
M t will be performed. 

The first two-loop results for Ap in the limits Mh — > and Mh 3> m t for large m f can 
be found in j2], for arbitrary Mh in [3]. The result for Ap of order G F m 2 a 2 was obtained 
in jlj. By including terms of C(M^/m 2 ) this has even been extended to predictions for 
Ar in three-loop approximation. 

The present paper is concerned with mixed QCD and electroweak contributions of order 
(G F m 2 ) 2 a s and purely electroweak contributions of order (G^m 2 ) 3 . It extends the result of 
[5] where the special case M H = has been considered. Using the technique of heavy mass 
expansion for the limit Mjj 3> m t and considering sufficiently many terms leads to stable 
predictions down to Higgs Mass values as low as twice m t . Ap is furthermore evaluated 
for Mh around m t , and with the help of expansions, a stable approximation is provided 
in the region from close to zero up to twice m t . Combining the results from the different 
regions, the three-loop corrections to the p parameter are well under control. 

An outline of the calculation, including a short discussion of the renormalization, is 
given in section |21 In section 0] the transformation of the top quark mass from the MS 
scheme to the on-shell top mass is treated in detail. Our results 1 at order G F and a s G F 
are presented in sections El and 03 in the MS and on-shell definition for the top quark mass 
respectively. 

2 Treatment of diagrams and renormalization 

The computation of the three-loop contributions to Ap at orders G 3 F and a s G F for Mh ~ 
m t and Mh 3> m t is similar to the one discussed in [Hj for Mh = 0. The same Feynman 
diagrams contribute: W and Z boson self energies up to three loops, top quark self energies 
up to two loops, one-loop Higgs boson self energies, and Higgs tadpoles up to two loops. 
All diagrams were generated using the program QGRAF [6 . 

In comparison to [3] a non-vanishing Higgs mass introduces a second mass scale into 
the problem. In order to calculate such diagrams we consider two different ranges of the 
Higgs mass: First, in the range of Mh around m t we perform a straightforward Taylor 

1 Only numerical results are presented for the values of the coefficients of the various expansions. Their 
representations in terms of fractions and transcendental numbers can be obtained from the authors. 
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expansion in the mass difference Mh — rn t around the point Mh = mt. In the second 
region, i.e. for Mh 3> fn t , we use asymptotic expansions |7j and made use of the program 
EXP jH]. In both cases the resulting diagrams contained at most one non- vanishing mass 
scale. 

From equation (J2J) the vector boson self energies can be taken at zero external mo- 
mentum. After the application of the expansion procedures the resulting integrals are of 
tadpole type with one mass scale. Their computation can be performed using the program 
package MATAD [9 j written in FORM [10J. Apparently the Higgs tadpoles, which are needed 
to renormalize the vacuum expectation value of the Higgs boson, can be treated using the 
same program package. To renormalize the mass of the Higgs boson, the computation of 
on-shell Higgs boson self energy diagrams is necessary. The corresponding renormalization 
constant is needed to one-loop order only and the relevant integrals are straightforward. 
The remaining two-loop top quark mass renormalization is, however, more involved, in 
particular in the on-shell scheme. Details on this issue are given in sections EH and EJ 

The whole calculation was performed using dimensional regularization and anticom- 
muting 75. This prescription preserves the Ward-Takahashi identities which relate the self 
energies of the gauge bosons to the non-diagonal self energies of the gauge-to-goldstone and 
the goldstone boson self energies. The validity of these identities was verified by explicit 
calculations. 



3 The p- parameter in the MS scheme 

In this section we present our results for the MS definition of the top quark mass. The 
renormalization constant of the top quark mass can be obtained by calculating the relevant 
one- and two-loop self energy diagrams in the limit of vanishing external momentum. After 
the application of the expansions mentioned in section El the remaining integrals were 
computed using the program package MATAD. 

The radiative corrections to Ap depend on the MS top mass m t = m t (m t ), the on- 
shell Higgs mass M H , the Fermi constant Gp and the strong coupling constant a s and are 
written as follows: 



Ap = x t Ap<**> + x 2 , Ap^?) + — x t Ap( QsX ') + ( —) 2 x t Ap 

7T ' V 7T / 



(a 2 s x t ) 



with 



+ ^ x\ Ap^ + x\ ApW + • • • , (3) 

7T 



= 3.0 • 10- 



Here and later on we always set the MS renormalization scale p equal to m t 

In this paper we evaluate the coefficients Ap^ x *> and Ap( asXt \ Numerically they are 
given by: 
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M H = 0: 



„3 



ApW = 47.73 , (4) 
Ap (a sX f) = 76 91 _ ( 5 ) 

Mh ~ mj: In this case we have computed five terms in the expansion in 5 defined by 
the relation Mh = m t (l + 5): 

Ap^ = -40.30- 138.26 5 -32.39 5 2 + 10.32 5 3 -3.275 4 - 0.01 5 5 + 0(5 6 ) , (6) 
Ap M) = 93.07 + 21.275 + 0.535 2 -2.295 3 + 1.145 4 -0.515 5 + O(5 6 ). (7) 

Mh 3> Tn t : The application of the hard mass procedure leads to an expansion in the 
ratio y = Am 2 /M H : 

ApW = -(-12.80 + 24.75 logy) 

y 

- 192.86 - 394.05 log y - 275.02 log 2 y - 24.25 log 3 y 

+ y (-337.88 + 408.59 log y - 233.35 log 2 y + 108.04 log 3 y) 

+ y 2 (82.56 - 378.31 logy + 26.14 log 2 ?/ + 48.60 log 3 ?/) 

+ y 3 (—82.58 — 291.03 logy — 11.37 log 2 y + 56.87 log 3 ?/) (8) 

+ y 4 (-44.06 - 366.54 log y - 3.27 log 2 y + 59.62 log 3 y) 

+ y 5 (-42.49 - 409.57 log?/ - 4.44 log 2 ?/ + 64.70 log 3 ?/) 

+ 0(y 6 ), 

A p M) = 216.07+ 131.38 log?/ + 57.07 log 2 ?/ + 9.00 log 3 ?/ 

+ y (-59.99 - 66.29 log y + 18.39 log 2 y - 19.22 log 3 y) 

+ y 2 (-66.31 + 21. 04 log?/ - 37.76 log 2 ?/ -5.57 log 3 ?/) 

+ y 3 (-9.73 - 28.56 log?/ - 29.38 log 2 y - 5.46 log 3 y) (9) 

+ y 4 (-4.48 -37.23 log?/ -30.13 log 2 ?/ -5.25 log 3 ?/) 

+ y 5 (-2.49-41.12 log?/ -31.92 log 2 ?/ -5. 14 log 3 ?/) 

+ 0(y 6 ). 

In Fig. [T]and El these results are plotted as functions of the ratio of the Higgs mass Mh and 
the top mass m t . The left scale gives the coefficient of the contribution under consideration. 
The right scale includes the prefactor for m t = 165 GeV, so the full numerical effect on Ap 
can be directly read off. In Fig. El the numerical value a s (m t ) = 0.109 was used for the 
right scale. 

The solid lines represent the sum of all terms displayed above. The dashed curves 
represent successive approximations. From Figs. Q and El we infer that the series given 
in equations (jUJ) to © provides a satisfactory approximation to the full result. In fact, 
a fairly smooth transition between the two approximations is observed for Mh ~ 2.5mt- 
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Figure 1: Contributions of order to Ap in the MS definition of the top quark mass. The 
black squares indicate the points where the exact result is known. 
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Figure 2: Contributions of order a s x\ to Ap in the MS definition of the top quark mass. 
The black squares indicate the points where the exact result is known. 
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In addition the expansion seems to provide a good approximation even down to Mh = 0, 
where it can be compared to the exact result. 

Note, that the dominant term of the three-loop electroweak corrections Ap( x *' in the 
limit of large Higgs masses is proportional to M H . This is in agreement with the screening 
theorem JIT] stating that the highest power correction depending on the mass of the Higgs 
boson is at least one power of Mfj below the expectation based on naive power counting. 
The terms ~ M H stem directly from the large mass expansion of the three-loop W and Z 
self energies. However, the absolute size is small when the prefactor is included. 

The smooth behavior of the two expansions in Figs. QandEltogether with the agreement 
with the result for Mh = leads to the assumption that the three-loop corrections are well 
under control in the MS scheme. 



4 Transition from MS to on-shell top quark mass 

In a final step the result for the on-shell definition of the top quark mass are presented. 
This requires the ratio between the MS renormalized top quark mass m t and the on-shell 
mass M t in two-loop approximation, again in the limit of large top quark mass. Let us 
parameterize the ratio between MS in the following way 

m* = i + ^ cf(«.) + Xt C^ + C&> + ^X t C^ + X*Cl*n + . . • , (10) 

M t TT \ TT J TT 

using the expansion parameter 

X i = ^S^3.2-10- 3 , 
8V2tt 2 

which depends on the on-shell top quark mass M t . 

The results for (e.g. |T2! and [T3] respectively), and [Uj are well 

known. Results for both C^ asXi ^ and C^ x *' were obtained from the large mass and the 
mass difference expansion and are needed for the present discussion. 

The corresponding two-loop top quark self energy diagrams, already discussed in sec- 
tion El for vanishing external momentum, have to be recalculated for external momentum 
on mass shell. In general the relevant diagrams can be divided into two classes: those 
containing at least one virtual Higgs boson and those without Higgs bosons. The latter 
have already been evaluated for the computation of Ap in the case Mh = jSJ. However, 
the former class of diagrams needs additional attention. 

Let us first consider the limit Mh ^> m t . Asymptotic expansions reduce the two- loop 
on-shell integrals containing two mass scales to combinations of one scale integrals. These 
one scale integrals include: tadpole integrals up to two loops, one-loop on-shell integrals, 
and one and two-loop integrals containing a small (compared to the internal Higgs mass) 
external momentum. Each of these integrals is either straightforward to calculate or can 
be treated using MATAD and MINCER [15 . Thus one obtains the on-shell top quark self mass 
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in terms of the MS mass in two-loop approximation for the region Mh 3> fn t . The result 
reads as follows: 

C (x ' 2) = -(1.61 - 18.00 log K) 

+ 8.09 — 2.10 logF — 1.50 log 2 Y 

+ Y (-4.50 + 6.05 log Y- 1.73 log 2 Y) 

+ Y 2 (14.60 - 9.74 log Y + 1.69 log 2 Y) (11) 

+ Y 3 (2.02 - 2.34 log Y + 0.64 log 2 Y) 

+ Y 4 (0.56 - 0.94 log Y + 0.32 log 2 Y) 

+ F 5 (0.22 — 0.48 logF + 0.17 log 2 Y) + 0(Y 6 ) , 

C (*.X t ) = 7 09 _ 3 gg 1()g y + L5Q bg 2 y 

+ F (-3.78 + 2.93 log Y - 0.36 log 2 Y) 

+ F 2 (—0.89 + 0.99 log Y — 0.23 log 2 Y) (12) 

+ Y 3 (-0.28 + 0.42 log Y- 0.15 log 2 Y) 

+ F 4 (-0.11 + 0.21 logy- 0.11 log 2 Y) 

+ y 5 (0.22 - 0.48 io g y + o.i7 log 2 y) + c(y 6 ), 

where y = AM 2 /M H is used. 

For M# around m t two- loop on-shell diagrams composed of several massive and massless 
lines enter. Let us start with the case Mh = m t which involves two-loop on-shell diagrams 
with internal propagators which are either massless or of mass m t . The program packages 
0NSHELL2 jTHj and TARCER [T7j . based on the recurrence relations [TSj . are specifically 
tailored to this case. In the next step the Taylor expansion around this point is performed 
and five terms of order a s x t and x\ are evaluated. It should be mentioned, that in order 
x\ there are several diagrams with threshold at M H . This means that in order to expand 
these diagrams the threshold expansion should be applied EJ] . By explicit calculation, 
however, we find that the first five coefficients are given by the naive Taylor series. The 
results read as follows: 

= 33.50 - 22.91 5 + 13.96 5 2 - 2.30 5 3 + 2.71 5 4 - 3.46 5 5 + 0(5 6 ) , (13) 
C (a 8 x t ) = _22.15 + 34.76 5 - 6.83 5 2 + 2.64 5 3 - 1.32 <5 4 + 0.19 5 5 + (5 6 ) , (14) 

using S defined by M H = M t (l + 5). 

For completeness we also include the result for Mh = 0: 

C^) = 83.58 and C {(XsXt) = -17.53. (15) 

The results are shown in Fig. El and EJ A smooth transition between the two regions 
of large Mh and Mh ~ M t is observed. Furthermore, the expansion around Mh ~ M t , 
when extrapolated down to Mh = 0, leads to a remarkably good agreement with the value 
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Figure 3: Contributions of order xj to the relation between the MS and on-shell top quark 
mass. The black squares indicate the points where the exact result is known. 
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obtained from an analytical calculation at Mh = 0, a further and independent test of the 
approximation procedure. 

The result for the ratio between MS mass and on-shell mass — as well as the analytical 
expressions for the coefficients given in sections El and — is available from the authors 
upon request. 



5 Ap expressed through the on-shell top quark mass 

Combining the results of sections El and 0] we compute Ap in terms of the on-shell top 
quark mass M t , the on-shell Higgs mass Mh and the low energy Fermi constant Gp. Our 
calculation reproduces the one- and two-loop results for the p parameter at orders Gp, G% 
and G F a s [HI2I3]. 

The result for the on-shell definition of the top quark mass is parameterized in the same 
way as the MS result in section 01 In order to distinguish the two schemes we use again 
the expansion parameter X t , depending on the on-shell top quark mass M t . 

For completeness we also include the results for Mh = 0, computed in 

Ap (x ' } = 249.74 , (16) 

For the expansion around Mh = M t we obtain: 

Ap^ = 95.92 - 111.985 + 8.0995 2 + 9.365 3 + 7.27<5 4 + -15.60 5 5 + 0(5 6 ), (18) 
Ap (aJ '» = 157.295 + 112.005 - 24.735 2 + 7.395 3 - 3.525 4 + 2.06 5 5 + 0(5 6 ), (19) 

using 5 defined by M H = M t (l + 5). 
In the limit Mh 2> M t we obtain 

Ap (x ^ = ^(-3.17- 83.25 logy) 

- 189.93 - 231.48 log Y - 142.06 log 2 Y + 2.75 log 3 Y 

+ Y (-332.34 + 77.71 log Y - 68.67 log 2 Y + 51.79 log 3 Y) 

+ Y 2 (227.55 - 510.55 log Y + 87.77 log 2 Y + 6.41 log 3 Y) 

+ Y 3 (-58.40 - 329.18 log Y + 20.42 log 2 Y + 14.54 log 3 Y) (20) 

+ Y 4 (-36.14- 381.88 logF + 18.63 log 2 Y+ 15.04 log 3 Y) 

+ F 5 (-39.08 -416.36 log Y + 13.76 log 2 Y + 17.19 log 3 Y) 

+ 0(Y 6 ), 

Ap Kx t 2) = 7g 73 _ 47 _ 77 log y + 42.07 log 2 Y + 9.00 log 3 Y 

+ Y (225.16 — 179.74 logF + 70.22 log 2 Y — 19.22 log 3 Y) 

+ Y 2 (-76.07 + 25.33 logF- 9.17 log 2 Y- 5.57 log 3 Y) (21) 
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Figure 5: Contributions of order Xf to Ap in the on-shell definition of the top quark mass. 
The black squares indicate the points where the exact result is known. 
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+ Y 3 (—10.10 — 24.69 logF — 0.30 log 2 Y — 5.46 log 3 Y) 

+ Y 4 (-4.52 - 32.85 log Y + 0.72 log 2 Y — 5.25 log 3 Y) 

+ Y 5 (—2.55 — 36.61 logF + 1.06 log 2 Y — 5.14 log 3 Y) 

+ 0(Y 6 ), 



with 

AM 2 

Y = ^-. (22) 

The results in the on-shell definition of the top quark mass are shown in Fig. 03 and 
IH1 Again, the solid line corresponds to the numerical values including all terms of our 
expansion. The left and right scales give the size of the coefficient and the full effect on 
Ap respectively. 

The a s Xf contributions are sizeable and monotonously increasing over the full range of 
the Higgs mass suggesting a smoothly increasing behavior qualitatively similar to the well- 
known X 2 correction (displayed in Fig. |7J). For Mjj = the result seems to be accidentally 
small, but increases rapidly for larger Higgs mass values. The Xf term exhibits a minimum 
in the range between Mh = 2M t and Mh = 3M t and appears to be very small or even 
negative in this region. 

To compare the size of the three-loop contributions to the known one- and two-loop 
terms we display their individual effects on Ap in Fig. including the corrections of 
order a 2 X t computed in [I]. We start with the sum of the contributions of order X t 
and a s X t which are independent of the Higgs mass and proceed by adding the particular 
orders to Ap. To identify the individual contributions for Mh = we continued the value 
horizontally to the negative Mn/M t region. 

Let us first consider the a s X 2 and Xf terms for the case Mh = 0, discussed already in 
Inclusion of the a s X 2 term leads to a change in Ap of 1.0 • 10~ 6 , which is not visible 
in Fig. [7| The Xf term increases Ap by roughly 8.2 ■ 10 -6 , again an extremely small 
effect. The picture alters drastically for larger values of Mh- The order a s X 2 correction 
is sizeable in the region of Mh ~ M t and above. It amounts to approximately 1/3 of the 
two-loop contribution Ap^ Xt \ but appears with the opposite sign. In contrast, the pure 
electroweak three-loop contribution Ap( Xt ) yields an enhancement of about 1% compared 
to the order X 2 terms and is therefore negligible for Mh values of present interest. These 
corrections to Ap can be used to examine the resulting shifts of observables like the W 
mass and the effective weak mixing angle, predicted for fixed values of M z , a, and Gp 
through the equations 




^-^i^-^G^-) (23) 



and 



Mjjr 1 /. L 4vra 




sta ^= 1 -^ = H 1 -V 1 -V5^lT^ + ---)l> (24) 
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where Aa is the shift of a due to photon vacuum polarization effects and the dots are 
nonleading remainder terms. In Fig. |S] the resulting shifts are shown for the different 
orders of Ap in dependence of the Higgs mass (We do not display the shift of order a s X t 
which amounts to —52.8 MeV). This shift should be compared to 5M\y = 6 MeV and 
5sin 2 # e fj = 1.3 ■ 10 -5 which is the precision anticipated for TESLA operating at the Z- 
resonance |2Hj . 

Alternatively we also present the shift in M t which would lead to a comparable shift in 

(25) 

M t 2 3X t 1 ; 

For 5Ap = 5 • 10~ 5 this leads to a shift of m t by about 500 MeV, which has to be compared 
to 5m t ~ (9(100MeV) anticipated for measurements at a future linear collider. 

For fixed M t , Mw, and Mz our result for Ap has to be compensated by a shift in the 
prediction for M#\ 

— — w — - . 2 5Ap w 670 5Ap 26 
% llGp sin v w 

For Mu between 120 GeV and 200 GeV the term a s X^ leads to an upward shift between 
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Figure 8: Shifts in My/ and in sin 2 6w from the various corrections to Ap. 



4.5 and 7.5 GeV, the Xf term to further positive shift of less than 0.5 GeV. 

The X^a s term is, furthermore, comparable or larger than the non logarithmically 
enhanced light fermion electroweak two-loop corrections j^I] and significantly larger than 
the purely bosonic two-loop contribution. The Xf is comparable in size to the bosonic 
two- loop effects evaluated in [22] . 

Summary: Combining the results from the extreme regions M# = 0, Mh ~ M t , and 
Mh 3> Mt a prediction for the Xf and a s Xf contributions to the p parameter has been 
obtained, which covers in a reliable way the full region of Mh- The magnitude of these 
terms is comparable to the non-enhanced two-loop corrections. 
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